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In their comment pi on our paper H], Ben-Naim and 
Krapivsky claims to have presented numerical evidence 
against our conclusions. We show that their claims are 
not valid for the following reasons, which we discuss here, 
(i) The arguments given in pi against non-universality 
are incorrect. 

(ii) The apparent disagreement between the results arises 
possibly from a difference in the initial conditions used, 
(iii) We present new numerical results, which support our 
earlier conclusions. 

In [|l|], the authors emphasise that their scaling func- 
tion is independent of all initial conditions. This is a 
surprising result, because unlike the case of particle den- 
sity n{t), the prefactor T appearing in the asymptotic 
expression P{t) = Tt~^ depends on the initial density 
uq. To show this, we first note that for an initial particle 
density no, the effects of the annihilation reaction will be 
felt only after a time interval Iq ~ (DnQ)~^. For t ^ to, 
n{t) - t~i/2 and so P{t) = P{tQ){t/to)-^ . On the other 
hand, for t < to, n{t) « tiq so that P{t) ~ p(o)e-""»^^^ 
in this regime (following an argument given in |0]). Here 
P(0) is the initial density of persistent sites and a is 
a numerical constant. Combining the two, one finds 
r ~ P(0)no~^^, which we have verified in simulations. 
(Similar arguments can be used to show that the prefac- 
tor for particle density n{t) is independent of no). 

It is easy to see that 1/2 cannot be a universal value 
for the dynamical exponent z proposed in f^, contrary 
to the claims in [Q. To show this, let us consider the 
q-state Potts model in the limit q —> oo, the dynamics 
of which is given by A + A ^ A reaction which also has 
n{t) ^ t^^/^. It is known from the exact solution y] 
that 6 — 1 for this model. If we now follow the argu- 
ments in M, we find that their scaling function !F{x) ^ 
constant for x <C 1 and decays exponentially for a; ^ 1 
(where x = lt~'^''^). Recomputing the persistent fraction 
P{t) using this distribution yields the inconsistent result 

Pit) - i-l/2. 

Coming to the numerics, first of all, we would like to 
point out that L{t) curves in m for no — 0.5 and no = 0.8 
are almost identical, unlike the corresponding plots given 
in our paper. This is possibly due to a difference in the 
initial condition: the entire lattice being taken as persis- 
tent at t = i.e, P{0) = 1 instead of -P(O) = 1 — no as in 
M. In the former case, the fragmentation of the lattice 
into clusters of persistent and non-persistent sites is de- 
layed by one MC step (for D = 1/2). Consequently, the 
effective particle density is only no* = no — Sno, where 
dno is the decrease in density over one time step. We 
have found numerically that for no = 0.2, 0.5 and 0.8, 
no* — 0.15, 0.25 and 0.29 respectively. In this small 


range, it is difficult to make any inference regarding the 
universality of exponents. 
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FIG. 1. L{t) is plotted against t for no — 0.95 and 
no = 0.8. For comparison, the plot for no = 0.8 with the 
initial condition P(0) = 1 is also shown, which is identical to 
the no = 0.8 plot in the comment. Our system size is 10^ 
and the results were averaged over 100 different initial real- 
izations. We have also done one 10^ size lattice, with identical 
results. Note that for large no, the asymptotic behaviour sets 
in very early. 


Finally, we present more numerical results (Fig. 1) 
which endorses the conclusions in our paper. In particu- 
lar, we see that as no — > 1, z c 
with our earlier observations 


— (Table I), in agreement 


no 
0.8 
0.95 


0.37515(1) 
0.37573(1) 


0.44916(2) 
0.39517(2) 


0.46894 
0.39551 


TABLE I. Exponents measured from simulations on 10'' 
size systems. Figures in brackets represent statistical error in 
the last decimal place. 
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